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Abstract 

Diffusion of electrons in two-dimensional disordered systems with spin-orbit 

interactions is investigated numerically. Asymptotic behaviors of the second 

moment of the wave packet and of the temporal auto-correlation function are 

examined. At the critical point, the auto-correlation function exhibits the 

power-law decay with a non-conventional exponent a which is related to the 

fractal structure in the energy spectrum and in the wave functions. In the 

metallic regime, the present results imply that transport properties can be 

described by the diffusion equation for normal metals. 
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The Anderson transition in disordered systems is usually classified into three universality 
classes, the orthogonal, the unitary and the symplectic [|I|]. The classification depends only 
on the symmetry of the system. In the case where the Hamiltonian is represented by a 
real and symmetric matrix, the system belongs to the orthogonal class. Systems having 
no time reversal symmetry belong to the unitary class. The system with a spin-dependent 
interaction having the time reversal invariance belongs to the symplectic class. 

In two dimensions, these three universality classes show quite different critical behaviors. 
In the case of the orthogonal class, all the states are considered to be localized. In the case of 
the system in a magnetic field, which is a typical example of the unitary class, all the states 
are again localized except at the centers of the Landau bands. Anyway, in two dimensions, 
these two classes do not show the Anderson transition. In contrast, the symplectic class is 
known to exhibit the Anderson transition in two-dimensional systems 

A typical system which belongs to the symplectic universality class is a disordered system 
with spin-orbit interactions. It has been proved that in two dimensions, such a system 
has localized states if the disorder is strong enough compared to the spin-orbit coupling 
Numerical and analytical studies including the finite-size scaling method and the 
renormalization group theory P|-pT| have been made on the critical behavior of the Anderson 



transition in two-dimensions. However, several important issues still remain to be clarified. 
In particular, the transport properties in the metallic regime appear to be peculiar. First, 
a direct application of the scaling theory leads to a result that the conductivity diverges 
logarithmically |p|,pT|. The infinite conductivity in the presence of the scattering by the 
impurities does not seem to be physically acceptable. Second, the scaling function in the 
metallic regime, obtained by the finite-size scaling method MjM, is different from what we 



expect from the theory developed by MacKinnon and Kramer |T2[ for the orthogonal case. 
These facts imply that the metal-insulator transition in two dimensions may have some 
complex structures which have to be explored further. The validity of the one-parameter 



scaling hypothesis in such a case, for instance, has already been discussed [QII^- It is thus 
an interesting problem to clarify the transport properties in the metallic regime and at the 



critical point of two-dimensional systems. 

In this letter, we examine the asymptotic behaviors of the temporal auto-correlation 
function and of the second moment of the wave packets in the two-dimensional disordered 
system with spin-orbit interactions. Since the Anderson transition in two dimensions occurs 
only in the case of symplectic ensemble the emphasis is placed on the behaviors of these 
quantities in the metallic regime and at the critical point. The second moment of the wave 
packets and the temporal auto-correlation function are defined by 



{r\t))= J drr'\^p{r,t)\' (1) 

and by 

c{t) = i fp{t')dt', p{t) = \{m\m)\^ (2) 

t Jo 

respectively. Here ilj{t) and P{t) denote the wave function at time t and the overlap function 
between the initial state \ip{0)) and the state \ip(t)) at time t, respectively. 

In the metallic regime, the second moment in d dimensions is usually expected to grow 
in proportion to time t as 

(r^) ~ 2dDt, (3) 

where the diffusion coefficient is denoted by D. The diffusion coefficient D is related to the 
conductivity a by the Einstein relation a = e^Dp ||14[, where p is the density of states at 
the Fermi energy. 

The auto-correlation function C{t), on the other hand, is expected to decay in proportion 
to t^" in the metallic regime and at the critical point. It has been shown that the exponent 



a is related to the correlation dimension D2 of the energy spectrum as a = D2 p5| , p!6 
It has been further demonstrated that the correlation dimension D2 is also related to the 
generalized dimensionality D2 of the eigenfunctions as D2 = D2/2 p6|JT7[] . We can thus 
obtain information about the fractal structure in energy spectrum and in the eigenfunctions 
by evaluating the exponent a. 



If the time-evolution of the wave packet can be described by the diffusion equation, the 
auto-correlation function would decay in proportion to where d denotes the dimension- 
ality of the system. Thus the conventional value of a is 1 in two dimensions, and this value 
of a implies that D2 = 1 and D2 = d = 2, indicating the absence of any fractal structure in 
this system. The deviation of a from 1 can therefore be regarded as evidence for a fractal 
structure in energy spectrum and in the eigenfunctions. 

In the localized regime, the wave packet does not diffuse. Hence, both the second moment 
and the auto-correlation function have a finite and nonzero value in the limit t — ^ 00. 

The model we adopt in the present letter is the Ando model described by the Hamil- 
tonian 

H = J2e.Cla^- E ^-.vCtQ.' (4) 

i,cr <i,j>,a,cr' 



with 



y -V2 J 



^ Vi -iV2 ^ 



(5) 



where Cj^{Cif^) denotes a creation (annihilation) operator of an electron on the site i with 
spin 0" and Ei denotes the site-diagonal potential distributed independently. The probability 
distribution of the site-potential is uniform in the range [—W/2, W/2]. Here, the unit vector 
of the direction is denoted by x{y). All length-scales are measured in units of the 

lattice spacing a. The effect of the spin-orbit interaction is taken into account by the off- 
diagonal matrix element V2 which is chosen such that the relation V2/V = 0.5 holds, where 



V = yVi + V2- For the states at the center of the band {E = 0), the critical point of this 
model has been estimated as Wc = 5.74V by the finite-size scaling method PJ^. We confine 
ourselves in this letter to the case of the band center {E = 0) for simplicity. 

To solve the time-dependent Schrodinger equation numerically, we have used the method 
based on the exponential product formulas |I^-pO[]. One of the advantages of our method 



is that we can treat much larger systems than those treated in the exact diagonalization 
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method. It should be emphasized that in the present approach, the time-reversal sym- 
metry, which characterizes the symplectic ensemble, is exactly preserved if the symmetric 
decomposition formula is adopted. 



The essence of the present method is to evaluate the time-evolution operator U{t) = 
exp{—iHt/h) by using the decomposition formula for exponential operators [|TB|. The n-th 
order decomposition Un satisfies the condition 

= f/„(5t) + 0(5r+^), (6) 

where 6t denotes the time step of the simulation. We have adopted the forth-order decom- 
position, as in the previous paper [^, given by 



[/4 = U2{-ipt/h)U2i-i{l - 2p)t/h)U2{-ipt/h) (7) 



with 



where H = Hi + - ■ ■ + Hg and p = {2— v^)~^. The decomposition f/„ is called symmetric if it 



satisfies the condition that Un{t)Un{—t) = 1 ||T8[. It is easy to verify that the decomposition 



U4 indeed satisfies this condition. It can be shown that the symmetric decomposition U4 
has the property 

QU^itye-^ = Uiit), (8) 

provided that the Hamiltonians Hi, . . . , Hq are time reversal invariant, where denotes the 
time-reversal operator Consequently, we have relationships for the time-evolution of 

the states as {f, +\U4{t)\j; -) = -{j; +\Ui{t)\i; ~) and {i; +\Ui{t)\j; +) = {j; -\Ui{t)\i; -). 
Here \i; ±) represent the states of the spin up(+) and down(— ) electrons at the site i. It 
should be emphasized that preserving the symmetry of the system in the numerical approach 
is quite important because the properties of two-dimensional disordered systems strongly 
depend on their symmetries. 
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Numerical calculations are performed in systems whose sizes are 399 x 399 ~ 999 x 999 
and for 5 ~ 10 realizations of potentials. The initial state is built by diagonalizing a system 
of radius L = 21 located at the center of the whole system I^O],^. The time step is chosen 
as 6t = 0.2n/Vi. 

Let us first discuss the temporal auto-correlation function. The logarithm of the auto- 
correlation functions logiQ C{t) for the localized regime {W = 9V), the extended regime 
{W = IV) and the critical point {W = 5.7AV) are shown in Fig. 1. For W = IV, we find 
that the auto-correlation function exhibits the power-law decay and its exponent is evaluated 
to be 1. This implies that the system does not have any fractal structure here and that the 
diffusion of electrons in this region can be described by the diffusion equation. At the critical 
point W = Wc, the auto-correlation function again shows the power-law decay. However, 
in contrast to the case of W = IV, the exponent takes a value different from 1, reflecting 
the fractal structure at the critical point |[T7|j23[| in energy spectrum and in eigenfunctions. 



In the localized regime {W = 9V), the auto-correlation function is likely to have a nonzero 
value in the limit t oo, as expected. 

At the critical point, the value of the exponent a is estimated to be « = 0.84 ± 0.03. By 
assuming the relation D2 = 2a ||TR|T7|, the generalized dimension D2 of the wave function is 



then obtained as D2 = 1.68 ± 0.06. Our estimate of the exponent a agrees with the value of 
the correlation dimension D2{= 0.83 ± 0.03) obtained by the exact diagonalization method 



T7| and is consistent with the relationship a = D2 

Furthermore, we have made simulations for other types of distribution of site energies. 



namely, the Gaussian and the binary potential distributions p|,p!3|. The exponents at the 



critical points are estimated to be a = 0.85 ± 0.02 for the Gaussian distribution and a = 
0.85 ±0.03 for the binary distribution. Together with the result for the uniform distribution, 
the present results suggest the universal fractal dimensionality at the critical points of two- 
dimensional symplectic systems. 

We have also performed the simulations in the metallic regime for W = 3,4 and 5V. 
Here the auto-correlation function keeps to show the power-law decay. The estimated values 
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of the exponent a are summarized in Table 1 . As the strength of disorder approaches to the 
critical value the exponent starts to deviate from 1. This crossover phenomenon might 
be interpreted as a finite-size effect. As demonstrated in ref. [^], a fractal structure can be 



observed at length-scales much smaller than the characteristic length of the system. The 
characteristic length of the metallic regime, namely the correlation length, becomes larger 
than the size of the system as the critical point is approached and diverges at the critical 
point [^. It is thus expected that the exponent a deviates from 1 accordingly near the 
critical point, reflecting the fractal character of the present system at length-scales much 
smaller than the characteristic length. 

Let us consider, next, the behavior of the second moment of the wave packets. Its growth 
in the metallic regime for W = 1, 3, 4 and 5V are shown in Fig. 2. Here, the second moment 
increases linearly in time and thus, within the present length-scales, the conductivity stays 
at a finite value. This linear increase of second moment is also observed at the critical point 
as shown in Fig. 3. 

From our estimates of the diffusion coefficient and the density of states per area p ~ 
0.25/(a^y), the conductivity at the critical point {W = Wc) can be evaluated as dc ~ 
1.5(e^//i), which is consistent with the estimation in ref. |[TB|. In the localized regime, we 
observe that the growth of the second moment is apparently suppressed (Fig. 3). The values 
of the diffusion coefficient D estimated using the relation @ for various values of W in the 
metallic regime are also listed in Table 1. 

Fal'ko and Efetov have recently proposed a formula for the fractal dimensionality D2 in 
terms of the diffusion coefficient D in the weakly localized regime {27rpDh ^ 1) |2^. In 



the case of the symplectic ensemble, it is expressed as D2 = 2 — l/(27r^pD). Together with 
a = D2/2, our results (Table 1) are in good agreement with this formula. 

In summary, we have performed numerical simulation for the diffusion of electrons in 
the disordered two-dimensional system with the spin-orbit interaction. The present method 
preserves the time-reversal symmetry which characterizes the symplectic universality class. 
In the metallic regime away from the critical point, the asymptotic behaviors of the temporal 



auto-correlation function and of the second moment of the wave packets are nothing but 
those expected from the diffusion equation. Within our simulation, there is no evidence for 
the logarithmic divergence of the conductivity, which is a consequence of a direct application 
of the scaling theory P,pl]]. At the critical point, the temporal auto-correlation function 
shows the power-law decay with a non- conventional exponent reflecting the fractal structure 
in energy spectrum and in the wave function, while the growth of the second moment seems 
to be diffusive. We note that this type of behavior has been observed in the critical states 
in the quantum Hall system |jl6|] as well as in the system with random magnetic fields 
||20| . It would therefore represent a common property at the critical point of the Anderson 
transition in two dimensional systems. The value of the exponent a at the critical point is 
consistent with the correlation dimension obtained in ref. [0 and the relation a = D2 
|T5| . We have further evaluated the exponent a at the critical points for several types of 
the potential distributions and found that those values coincide with each other within the 
present numerical accuracy. This fact indicates the universality of the exponent a at the 
critical points in two-dimensional symplectic systems, which would enable us to determine 
the critical points by using the equation of motion method. 
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FIGURES 

FIG. 1. The temporal auto-correlation function C{t) for the extended regime {W = IV), the 

localized regime (W = 9V) and at the critical point {W = Wc). The exponent a is estimated as 
0.84 ± 0.03 and 0.996 ± 0.014 for W = Wc and W = IV, respectively. 

FIG. 2. The growth of the second moment in the metallic regime for 1^ = 1, 3, 4 and 5V. Solid 
lines correspond to D = 17.2, 2.87, 1.99 and lA5a^V/h for W = 1,3, 4, and 5^, respectively. 

FIG. 3. The growth of the second moment at the critical point {W = Wc) and in the localized 
regime (W = 8, 9V). Solid line for W = Wc corresponds to D = 0.97a V/^. 
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TABLES 



WjV D/{a'^V/n) p/(l/(aV)) a 

1 ~ 17.2 ~ 0.36 0.996 ± 0.014 

3 - 2.87 - 0.35 0.97 it 0.04 

4 ~ 1.99 ~ 0.30 0.96 ± 0.03 

5 ~ 1.45 ~ 0.27 0.93 ± 0.03 
5.74(= Wc/V) - 0.97 ~ 0.25 0.84 ± 0.03 



TABLE I. List of the exponent a, the diffusion coefficient D and the density of states p for the 
disorder W/V = 1,3,4,5 and 5.74(= Wc/V). Here a denotes the lattice spacing. 
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